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1. INTRODUCTION 

Multiferroic materials exhibiting ferroelectricity and ferromagnetism in the same phase have attracted 
much attention due to their basic physics and potential technological applications in recent years [1, 2]. 
Lithium tantalate (LiTaO3) and lithium niobate (LiNbO3) are well-known materials for ferroelectric, 
piezoelectric, acoustooptical, electro-optical and nonlinear optical applications [3, 4]. The interest in the use of 
piezoelectric materials as the wave propagation medium lies in the propagation of appearance. These waves in 
this case, spread in a resilient part (or acoustic) and power; hence the name electroelastic or electroacoustic [5]. 
Various types of transducers such as bulk acoustic wave (BAW) (Bulk acoustic waves) transducers, shear wave 
transducers and interdigital transducers IDTs (Inter-Digitized transducers) are reported for generation and 
reception of acoustic waves in SAW (Surface acoustic waves) devices. IDTs are widely used transducer in 
SAW devices and they are metallic comb shaped electrodes fabricated over a piezoelectric substrate [6]. Theory 
on generation and propagation of BAW in a SAW device with IDTs is well explained in [5, 7]. 
Some materials, when deformed, become electrically polarized.This effect is known as piezoelectricity [8, 9]. 
In this work, we searched for all Acoustic velocity and Attenuation coefficients on the level of Piezoelectric 
substrate (LiNbO3) and then we compared it with the results of Piezoelectric substrate (LiTaO3) ,to conclude 
wich of them is better in utilization for generating bulk acoustic waves. 
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Figure 1. LiINbO3 or LiTaO3 Crystals excited by transducer. 


2. PHENOMENOLOGICAL TENSORIAL PIEZOELECTRIC EQUATIONS 
The electrical induction and the electric field E in the piezoelectric materials are bound by: [8-11] 


D=eE+e'S (1) 
Where € is the permittivity tensor (F/m)? 
e.: piezoelectric tensor (c/m) 


S: strain tensor (Sij=1/2[0Ui/0Xj + 0Uj/0Xi]) 


The electric polarization of the medium under the effect of deformation, also involves the creation of stresses 
under the effect of an external electric field. This constraint is: 


T=C.S-e.E (2) 
where C : represents the elastic tensor (N/m?) 
In one form the stress tensor and the electric induction are defined as follows: [8-11] 

Ty, = Cina Sia — ny Ex (3) 


DD, =¢54-Sy +64. Ey (4) 


1 


with i,j,k,1 =1,2,3 
In the quasi-static approximation, the Maxwell equations reduce to the Poisson equation: 


_ D. 
‘vie 20 (5) 
ox 


1 
The movement of the particles under the action of stress is described by Newton's equation: 


ou? 


VT=09. 
ot? 


(6) 


Substituting (3) and (4) in (5) and (6) we obtain the piezoelectric tensor phenomenological equations: 
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3. GENERAL FORM OF THE SOLUTION 

The piezoelectric wave is the solution of the elastic and electric equations that satisfy elastic and 
electric boundary conditions. When these waves attenuate exponentially inside the piezoelectric material as 
shown in Figure 1, they are known as surface waves and if this is not the case they are known as bulk waves. 
Consider the following form of the surface wave (partial wave). 


We = u,;,©expGB.a;, <>) (9) 


where ui(i=1,2,3) are displacement amplitudes, ui(i=4) represents the amplitude of the electric 
potential, B is the propagation constant, ai are the attenuation coefficients of the wave within the piezoelectric 
crystal (axis Y show in Figure 1) and @ is the angular pulsation. We will be interested in the coefficient ai. With 
a chosen LiNbO3 (Lithium Niobate) or LiTaO3 (Lithium Tantalate). Solutions of this type of wave correspond 
to waves propagating with or without attenuation along the direction X. The elastic displacements Ui and the 
electric potential U4 may vary from the normal direction to the flat surface (Y), but following invariant Z as 
show in Figure |. Substituting (11) into (9) and (10) we obtain the system: 


[A}[U] 7 [9] (10) 


Where [A] is a symmetric matrix 4x4 


[U] = [ul, u2, u3, u4] E Components determine 


3.1. For LiNbO3 (Lithium Niobate) 
Development of the determinant of [A] of the general, an eighth order polynomial: 


8 F 
Det(A) = 2 B.a = 0 (11) 


This polynomial is called dispersion equation or secular equation. 
The solution of (11) gives for each B (6 = 2.2.f/ Vs: constant propagation) eight roots. These roots are based 
on Vs (acoustic velocity). Each root generates three displacement components of the particle and an electric 
potential. Thus, the general solution is a combination of eight roots (8 secondary wave) given by the expression: 


8 
U, = S.C,.D?.exp{jp.a”.X, + jBX, (1+ 5-7) - 5B.V.t} (12) 
n=l 


i=1,4 
where Dare the components of the eigenvector of the system (10) associated with the eigenvalue q" 
C,: constant to be determined by the boundary conditions. 


3.2. For LiTaO3 (Lithium Tantalate) 
Development of the determinant of [A] of the general, an sixth order polynomial: 
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This polynomial is called dispersion equation or secular equation. 
The solution of (13) gives for each B (8 = 2.2.f/Vs: constant propagation) six roots. These roots are based on 
Vs (acoustic velocity). Each root generates three displacement components of the particle and an electric 
potential. Thus, the general solution is a combination of six roots (6 secondary wave) given by the expression: 


8 
U, = S.C, .D?-exp{jB.a".X, + jBX,(1+ jy) - i.V.t} (14) 
n=l 


i=1,3 
where Dare the components of the eigenvector of the system (10) associated with the eigenvalue " 
, Cn: constant to be determined by the boundary conditions. 


4. BEHAVIOR OF THE ROOTS OF THE SECULAR EQUATION 
The behavior of the roots depends on the acoustic velocity Vs at a velocity less transverse to the low 
velocity of volume in the medium, all the roots are complex conjugate pairs: 


im 


where q) = q{*) et a? =a"? with j=1,3,5,7 


This corresponds to the ‘partial wave: 


U® = exp j. B.0% .X, -p.a.x si exp(j6.X, ) (15) 


U® = exp j.p.a2.X, +p.a0.X , bexp(if.X,) (16) 
After development we have: 


U0 = exp|- —p.a0 .X | exp j.lal.X, + X, | (17) 


UO? = exp) +B. a xX = exp j.la.X, + X,] (18) 


Remarks 
if @=0 and at") <0, the bulk waves will be obtained [12-14]. 


5. PROBABILISTIC NEURAL NETWORK 

A Probabilistic Neural Network (PNN) is defined as an implementation of statistical algorithm called 
Kernel discriminate analysis in which the operations are organized into multilayer feed forward network with 
four layers: input layer, pattern layer, summation layer and output layer. A PNN is predominantly a classifier 
since it can map any input pattern to a number of classifications. [15, 16] Among the main advantages that 
discriminate PNN is: Fast training process, an inherently parallel structure, guaranteed to converge to an 
optimal classifier as the size of the representative training set increases and training samples can be added or 
removed without extensive retraining. Accordingly, a PNN learns more quickly than many neural networks 
model and have had success on a variety of applications. Based on these facts and advantages, 
PNN can be viewed as a supervised neural network that is capable of using it in system classification andpattern 
recognition [16]. The PNN consists of nodes allocated in three layers after the inputs show in 
Figure 2. [15, 17] 
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Summation Layer 


Output Layer 


Figure 2. The architecture of PNN [15,17] 


6. APPLICATION ON LiNbO3 (Lithium Niobate) 

In this example, we will show the variations of the real and the imaginary parts of attenuation 
coefficients (a3) show in Figure 3 and (a4) show in Figure 4, we will indicate the presence of bulk waves (zero 
imaginary part and real part is negative) and for a good detection of bulk waves we should have (Imaginary 


part equals zero and the minimal negative value of real part). We set the value of the real part for (a3), which 
we will look for it in the Classification as shown Figure 3. 


Real part 


zero imaginary 
part and the 
minimal negative a 


value of real part 


-1,1 -0,8 0,5 -0,2 0,1 0,4 
imaginary part 


Figure 3. Plot - Real part versus Imaginary part 


We set the value of the real part for (a4), which we will look for it in the Classification show in Figure 4. 
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Figure 4. Plot - Real part versus Imaginary part 


7. RESULTS AND DISCUSSION (for LiNbO3 Lithium Niobate) 
7.1. For coefficient attenuation a3 

To know the accurate values of each of the real part and the imaginary part which represent coefficient 
attenuation a3 and Acoustic velocity, Selection variable stays always Acoustic velocity and we change the 
Classification factor where we find: show in Figure 5, we have classification where we find classification factor 
is imaginary part and selection variable is acoustic velocity and number of cases in training set 51, in the 
minimal negative value of real part we find that the imaginary part equals 0. Show in Figure 6, we have 
classification where we find classification factor is real part and selection variable is acoustic velocity and 
number of cases in training set: 51, We find that the minimal negative value of real part equals -0.5654. Show 
in Figure 7, we have classification where we find classification factor acoustic velocity and selection variable 
acoustic velocity and number of cases in training set: 51, in the minimal negative value of part part (-0.5654) 
We find that the acoustic velocity equals 4050 m/s. 
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Figure 5. Neural Network Bayesian Classifier - Imaginary part (Acoustic velocity) 
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Figure 6. Neural Network Bayesian Classifier - Real part (Acoustic velocity) 
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Figure 7. Neural Network Bayesian Classifier - Acoustic velocity (Acoustic velocity) 


7.2. For coefficient attenuation a4 

To know the accurate values of each of the real part and the imaginary part which represent coefficient 
attenuation a4 and Acoustic velocity, Selection variable stays always Acoustic velocity and we change the 
Classification factor where we find: show in Figure 8, we have classification where we find classification factor 
is imaginary part and selection variable is acoustic velocity and number of cases in training set 51, in the 
minimal negative value of real part we find that the imaginary part equals 0. Show in Figure 9, we have 
classification where we find classification factor is real part and selection variable is acoustic velocity and 
number of cases in training set: 51, We find that the minimal negative value of real part equals - 0,57212. Show 
in Figure 10, we have classification where we find classification factor acoustic velocity and selection variable 
acoustic velocity and number of cases in training set: 51, in the minimal negative value of part part (- 0,57212) 
We find that the acoustic velocity equals 4000 m/s. 
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Figure 8. Neural Network Bayesian Classifier - imaginary part (Acoustic velocity) 
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Figure 9. Neural Network Bayesian Classifier - Real part (Acoustic velocity) 
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Figure 10. Neural Network Bayesian Classifier - Acoustic velocity (Acoustic velocity) 
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Table 1. Summarizes the full work where we can notice bulk waves detection. 
Coefficient Attenuation Imaginary part _ Real part __ Acoustic velocity _ Bulk waves (BAW) 


on All values All values All values No Detection 
Qo All values All values All values No Detection 
on 0 -0,5654 4050 m/s Good Detection 
O4 0 -0,57212 4000 m/s Good Detection 
Os All values All values All values No Detection 
O6 All values All values All values No Detection 
Qty All values All values All values No Detection 
Og All values All values All values No Detection 


8. APPLICATION ON LiTaO3 (Lithium Tantalate) 

In this example, we will show the variations of the real and the imaginary parts of attenuation 
coefficients (a2) show in Figure 11, we will indicate the presence of bulk waves (zero imaginary part and real 
part is negative) and for a good detection of bulk waves we should have (Imaginary part equals zero and the 
minimal negative value of real part). We set the value of the real part for (02), which we will look for it in the 
Classification show in Figure 11. 
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Figure 11. Plot - Real part versus Imaginary part 


9. RESULTS AND DISCUSSION (for LiTaO3 Lithium Tantalate) 
9.1. For Coefficient Attenuation a, 

To know the accurate values of each of the real part and the imaginary part which represent coefficient 
attenuation a3 and Acoustic velocity, Selection variable stays always Acoustic velocity and we change the 
Classification factor where we find: show in Figure 12, we have classification where we find classification 
factor is imaginary part and selection variable is acoustic velocity and number of cases in training set 51, in the 
minimal negative value of real part we find that the imaginary part equals 0. Show in Figure 13, we have 
classification where we find classification factor is real part and selection variable is acoustic velocity and 
number of cases in training set: 51, We find that the minimal negative value of real part equals -0.5654. Show 
in Figure 14, we have classification where we find classification factor acoustic velocity and selection variable 
acoustic velocity and number of cases in training set: 51, in the minimal negative value of part part (-0.5654) 
We find that the acoustic velocity equals 4050 m/s. 
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Figure 14. Neural Network Bayesian Classifier - Acoustic velocity (Acoustic velocity) 
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Table 2. Summarizes the full work where we can notice bulk waves detection. 
Coefficient Attenuation Imaginary part __ Real part Acoustic velocity _ Bulk waves (BAW) 


on All values All values All values No Detection 
Oy 0 -4,594 3450 m/s Good Detection 
fon All values All values All values No Detection 
4 All values All values All values No Detection 
Os All values All values All values No Detection 
O6 All values All values All values No Detection 


CONCLUSION 
In this article, we explained that the phenomenon of bulk waves are relying on numerical results at 


the level of attenuation coefficients. Changes in real and imaginary parts of the coefficients based on the 
acoustic velocity to detect these waves. Since the bulk waves were detected in Lithium Tantalate (LiTaO3) in 
one coefficient attenuation a, and in Lithium Niobate (LiNbO3) were detected in two coefficients attenuation 
a3 and a4 so we conclude that (LiNbO3) is the better to use in generating bulk acoustic waves. To compare 
between many piezoelectric substrates in detecting acoustic micro waves using PNN classification you just 
have to change the characteristics of the material parameters, and those results we need them in realization of 
acoustic microwaves devices. 
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